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STRONG SOLUTIONS FOR TWO-PHASE FREE BOUNDARY PROBLEMS FOR 
A CLASS OF NON-NEWTONIAN FLUIDS 

MATTHIAS HIEBER AND HIROKAZU SAITO 


Abstract. Consider the two-phase free boundary problem subject to surface tension and gravita¬ 
tional forces for a class of non-Newtonian fluids with stress tensors Ti of the form Ti = —irl + 
^i(\D{v)\^)D{v) for i = 1,2, respectively, and where the viscosity functions satisfy /ri(s) G 

C®([0, oo)) and /ri(0) > 0 for i = 1,2. It is shown that for given T > 0 this problem admits 
a unique, strong solution on (0, T) provided the initial data are sufficiently small in their natural 
norms. 


1. Introduction and main result 

The free boundary problem for two-phase flows for Newtonian fluids with or without surface tension 
is nowadays rather well understood. We refer in particular to the articles [22], [23], [21], [20], [24], [1], 
[4], [3] and [26] describing the present state of research for the situation of sharp interfaces. 

In order to describe the problem in more detail, let A > 2 and Fq C K.^ be a surface which separates 
a region ni(0) filled with a viscous, incompressible fluid from Il 2 ( 0 ), the complement of ni(0) in 
The region 112(0) is also occupied with a second incompressible, viscous fluid and it is assumed that 
the two fluids are immiscible. Denoting by r(t) the position of Fq at time t, T{t) is then the interface 
separating the two fluids occupying the regions Di(t) and 
An incompressible fluid is subject to the set of equations 

p{dtu + u ■ Vu) = div T, 
div u = 0, 

where p denotes the density of the fluid and the stress tensor T can be decomposed asT = r—pl, where 
p denotes the pressure p and r the tangential part of the stress tensor of the fluid. For a Newtonian 
fluid, r is given by r = 2p,D{u), where D{u) = [Vm + (Vu)^]/2 denotes the deformation tensor and p, 
the viscosity coefficient of the fluid. 

In this article we consider a class of non-Newtonian fluids, where r as above is replaced by 

r = 2p{\D{ur)D{u) 

for some function p satisfying 

(1.1) pGC^{[0,oo)) and ^(0) > 0. 

In the special case of power law fluids, one has 

p{\D{u)\^) = ,. + 

for some d > 1 and constants v,l3 >{). If d < 2, the fluid is then called a shear thinning fluid, if d > 2 
it is called a shear thickening fluid. Fluids of this type are special cases of so called Stokesian fluids, 
which were investigated analytically for fixed domains by Amann in [6] and [7]. 
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The motion of the two immiscible, incompressible and viscous fluids is then governed by the set of 
equations 

' p{dtv + V ■ Vv) = divT — p^a^N, in 


divi; = 0 in n(t), 

— iTrir] = (jHrnr on r(t), 

(1.2) < H = 0 on r(t), 

Vt = V ■ nr on T{t), 

ii|t=o = Vo in 

r|t=o = To, 


where n(t) = idi{t) U 02(t) and Cat = (0,..., 0,1)^. Here, the normal field on T{t), pointing from 
Oi(t) into is denoted by nr{t, ■). Moreover, Vr{t, •) and Hr{t, ■) denote the normal velocity and 

mean curvature of r(t), respectively. Furthermore, ja denotes the gravitational acceleration and a the 
coefficient of the surface tension. 

In this article we suppose that the stress tensor T is given by the generalized Newtonian type 
described above, that is, for given scalar functions pi,p2 ■ [0,oo) —>■ R, the stress tensor T is given by 

T = XQ.i{t)Ti{v, tt) + xn^{t)T2{,v, tt), T^{v, tt) = -ttI + 2pi{\D{v)\'^)D{v), i = 1,2, 

and where \D{u)\'^ = . The function xd denotes the indicator function of a set 

D C and the density p is defined by p := Xni{t)Pi + XQ.2{.t)P2 for the densities pi > 0 of the i-th 
fluid. The system is complemented by the initial fluid velocity iiq, the given initial height function ho 
and flo as well as Tq which are given by 

Ho = \ To, Fo = {(a;', xn) \ x' £ x^ = hoix')}. 

In addition, we denote the unit normal field on Fq by no- The quantity |/] = |/](a;,t) is the jump of 
the quantity /, which is defined on H(t), across the interface F(t) as 

lfj{x,t)= lim {f{x + enr,t) - f{x - enr,t)} for a; £ F(t). 

£—)- 0 + 

The problem then is to find functions v, tt and T solving equation (1.2). 

Well-posedness results for the above system (1.2) in the case of Newtonian fluids and in the special 
case of one-phase flows with or without surface tension were first obtained by Solonnikov [28], [29], 
[30], Shibata and Shimizu [25], [26]. The case of an ocean of infinite extend and which is bounded 
below by a solid surface and bounded above by a free surface was treated by Beale [9], Allain [5], Tani 
[32], Tani and Tanaka [33], Bae [8], and Denk, Geissert, Hieber, Saal and Sawada [13] and Gotz [16]. 

Besides the articles cited already above, the two-phase problem for Newtonian fluids was studied 
by Denisova in [11] and [12], and by Tanaka in [31] using Lagrangian coordinates. Indeed, Denisova 
proved local wellposedness in the Newtonian case in for r £ (5/2,3) for the case that one of 

the domains is bounded and Tanaka obtained wellposedness (including thermo-capillary convection) 
in for r £ (7/2,4). 

Priiss and Simonett were using in [21], [22] and [23] a different approach by transforming problem 
(1.2) to a problem on a fixed domain via the Hanzawa transform, which then was followed then by an 
optimal regularity approach for the linearized equations. Like this they proved wellposedness of the 
above problem in the case of Newtonian fluids. 

For an approach to the linearized problem based on Lagrangian coordinates , also in the setting of 
Newtonian fluids, we refer to the work of Shibata and Shimizu [27]. 

Problems of the above kind for non-Newtonian fluids were treated by Abels in [1] in the context of 
measure-valued varifold solutions. His result covers in particular the situation where p.i{s) = 
for 1 = 1,2 and d £ (1, oo). Note, however, that his approach does not give the uniqueness of a solution. 
For further results we refer also to the work of Abels, Dienig and Terasawa in [2]. Gotz investigated in 
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[16] the spin-coating process for generalized Newtonian fluids and showed local wellposedness of this 
problem for the setting of one-phase flows. 

Bothe and Priiss gave in [10] a local wellposedness result for Non-Newtonian fluids on fixed domains 
for viscosity functions /x e C^(0,oo) satisfying /i(s) > 0 and fi{s) + > 0 for s > 0. Note that 

our assumptions on the viscosity function ^ are different from those treated in [10]. For further results 
on various classes of non-Newtonian fluids on fixed domains we refer e.g. to the articles [14], [15], [17] 
and [19]. 

In our main result we show that system (1.2) admits a unique, strong solution on (0, T) for arbitrary 
r > 0 provided the viscosity functions Hi fulfill (1.1) and the initial data are sufficiently small in their 
natural norms. More precisely, we have the following result. 

Theorem 1.1. Let N + 2 < p < oo and J = (0, T) for some T > 0. Suppose that pi > 0, p 2 > 0, 7 a > 
0 , (7 > 0 and that 

Pi € C'^([0,oo)) and Pi{0) >0, i = 1,2. 

Then there exists eo > 0 such that for 

{vo^ho) € W'2-2/p(Oo)W X IF3-2/p(kN-1) 
satisfying the compatibility conditions 

lp{\D{vo)\'^)D{vo)no - {no ■ p{\D{vo)\'^)D{vQ)no}nQ] = 0 on Fq, 

div'(;o = 0 in LIq, |no] = 0 on Fq, 


as well as the smallness condition 


the system (1.2) admits a unique solution (n, tt, h) within the class 

V G H^piJ,Lp{nit))nLp{j,H^{n{t))f, 

7rGLp(J,i7i(0(t))), 

h G Lp(R^-i)) n VF2 -i/p(R^-1)) 

n IK V 2 -i/( 2 p) (1)) n Lp (J, w^-^^p (M^-^)). 

Remarks 1.2. a) Some remarks on notation are in order at this point. Setting 

= {{x\xn) 1x'gK^“\ a;A = 0}, 
by u G Lp{Q{t))) D Lp{J, Hp{Ll{t)))^ we mean that 

0*i; = u o 0 G H^{J, Lp(R^)) n Lp(J, 

where 0 and 0* are defined in Section 2 by (2.2) and (2.3), respectively. The regularity statement for 
TT is understood in the same way. 

b) The assumption p > TV -|- 2 implies that 

h G BUC{J,BUC^{R^-^)), dth G BUC{J, BUC^{R^-^)), 

which means that the condition on the free interface can be understood in the classical sense. 

c) Typical examples of viscosity functions p satisfying our conditions are given by 

p{s) = u(l -I- s^~ ) with d = 2,4, 6, or d > 8, 

d-2 

p{s) = v{l -I- s) 2 with 1 < d < oo 

for u > 0. For more information and details we refer e.g. to the work of [15], [14], [17] and [19]. 
Obviously, if d = 2, all viscosity functions corresponds to the Newtonian situation. 
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Let US remark at this point that our proof of Theorem 1.1 is inspired by the work by Priiss and 
Simonett in [22] and [23]. Our strategy may be described as follows: in Section 2 we transform the 
system (1.2) to a problem on a fixed domain. Maximal regularity properties of the associated linearized 
problem due to Priiss and Simoonett [23] are described in Section 3. Of special importance will be the 
function space F 3 (a) which will be introduced and investigated in Section 4. Finally, in Section 5, we 
treat the nonlinear problem and give a proof of our main theorem. 

In the following, the letter C denote a generic constant which value may change from line to line. 


2. Reduction to a fixed domain 


We start this section by calculating the divergence of the stress tensor, i.e. by calculating explicitly 

diY{fid{\D{u)\^)D{u)} for d=l,2. 

Let us remark first that given a vector u of length m for m > 2, we denote by Ui its i-th component 
and by u' its tangential component, i.e. u = (ui,..., and u' = {ui ,..., Um-i)'^■ We then obtain 

1 ^ 

{div{fid{\D{u)\'^)D{u)}), = 2 X! {‘^f^di\D{u)\‘^)D,j{u)Dki{u) + fJ,d{\D{u)\‘^)S,kSji}{djdkUi +djdiUk). 
For vectors u, v we set Ad{u)v := {Ad^i(u)v ,..., Ad,N{u)v)'^ where 


N 

Ad,iiu)v := - ^ A^J^/{D{u)){djdkVi + djdiVk) and, 

j,k,l=l 

:= ^(2fid{\D{u)\'^)Dij{u)Dki{u) + fid{\D{u)\'^)SikSji^ 

for d = 1, 2 and z = 1,..., iV. We then have 

Ad{u)u = — dW{fid{\D{u)\^)D{u)} and Al(i(0)u = —/id(0)(Au + V div m). 
In addition, we set 


A{u)v:= XQ^{t)Ai{u)v + Xn 2 {t)M{u)v and tt := tt + 
The system (1.2) may thus be rewritten as 

' p{dtv + V ■ Vv) — fi{0)Av + Vtt = —{A{v) — Al(0))z; 


( 2 . 1 ) 


div ?; = 0 

-|rnr] = aHrnr + IpjjaXN 

bl = 0 

Vt = V ■ nr 

v\t=o = Vo 

r|t=o = To; 


PlaXN- 


in fl(t), 
in fl(t), 
on F(t), 
on F(t), 
on F(t), 
in flo, 


where T = XQ.i{t)Ti{v,TT) + X^i 2 (t)^ 2 {v,TT) and p(0) = Xni(i)di(0) + Xn 2 (i)d 2 ( 0 ). 

Next, we transform the problem (2.1) to a problem on the fixed domain To this end, we define 
a transformation 0 on J x for J = (0, T) with T > 0 as 

(2.2) 0 : J X R^ 9 (r, ^ (L x', xm) € [J {s} x n(s), with t = r, x' = xn = , r) 


for some scalar-valued function h. Note that det^/0 = 1, where ^70 denotes the Jacobian matrix of 
0. We now define 


(2.3) 


u{x,0 ■= 0*v{t,x) := w(0(r,^)), 6>(r,0 := 0*7r(t,x) 
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as well as 

(2.4) 0,/(r, e) := /(0-'(x, t)) for , 

where 0“^ given by 0“^(a;,<) = {t,x',XN — h{t,x')). This change of coordinates implies the relations 

(2.5) dt = dr - {drh)DN, dj = Dj - {Djh)DN, djdk = DjDk - where 

J^jk(h) := {DjDkh)DM + {Djh)DMDk + {Dkh)DjDM — {Djh){Dkh)D^j^ 
for j,k = 1,... ,N, dr = d/dr and Dj = d/d^j since Dj^h = 0. 


Setting A' = V' = (9i,..., Sat-i)^, we first obtain similarly as in (2.5) 


(2.6) Dr{v) = E{u,h) := D^{u) — £{u,h), £{u,h) := (Dnu) 

Secondly, following [22, Section 2] we see that 

v^,h(t,e) 


'V^h 

T 

+ 

'V^h 

0 


0 


{DNuf. 


N-l 


i/r = E / ' ) = 

U Yi + iv',h(t,e)|2^ 
|V^,h|2A^,h 


= A^,h — 'H{h), where 


nih) := 


N-l 

E 


{Djh){Dkh){DjDkh) 


(l+^l + |V',hP)^l+|V',hP (l + |V',h|2)3/2 

Hence, system (2.1) is reduced to the following problem on 

pdrU — ^{Q)Au+ Vd = F{u,9,h) in R^, 
(i\\r u = Fd{u,h) in R^, 

-lp{0){DNUj + DjUn) 1 = Gj{u, I0j,h) on , 
- 2|/r(0)i:)ArUAr] - (Mya + aA')h = Gn^u, h) on R^, 


(2.7) 


|m] = 0 

drh — UN = Gh(u', h) 
w|t=0 = Uq 
h\t=o = ho 


on 


on 


on 


on 


j,A 

^0 > 


pAf 

^0 J 


t>N 


1 >N -1 


where j = 1,..., A — 1 and F = (Fi,..., Fn)'^■ The terms on the right hand side of (2.7) are given by 

N 

Fi{u, 0, h) := p{{drh)DNUi — (u • V)ui + {u ■ S/'h)DNUi} — p{0) + {Dih)DN0 + Ai{u, h) 

i=i 


= cr'H{h)Djh — {([p]7a + crA')h}Djh + {OjDjh + Bj{u, h) 
= —a'H{h)+BN{u,h), 

= [Dnu') ■ Vh = Dn{u' ■ V'h), 

= -u' • Vh. 


g,{uMM 

Gn{u, h) 

Fd{u, h) 

Gh{u, h) 

Here Ai{u, h), Bj{u, h) and Bn{u, h) are given by 

N 

A{u, h) := E h)) - ’^(0)) {DjDkUe + D^DiUk) 

N 


E [Al’'"’‘^iE{u,h))-AjE\0)YFjk{h)ue + Fje{h)uk), i 


= 1,...,A, 
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Bj{u, h) := - h)\^)DNUN'\Djh + l{fj,{\E{u, /i)p) - fi{0)}{DNUj + Djun)] 

N-l 

~ X! + DkUj)jDkh 

k^l 

N-l 

+ ^ lfi{\E{u, h)\‘^){DNUjDkh + DNUkDjh)jDkh, j = 1,..., - 1, 

k^l 

BN{u,h) :=2|{//(|E’('u,/i)p) -//(0)}Z)Ar'UAr] + lfj.{\E{u,h)\‘^)DNUNj\'^'h\‘^ 

N-l 

- ^ lfj,{\E{u, h)\‘^){DNUk + DkUN)jDkh 

k^l 

where 

Af'‘’\E{u,h)) := XR^N,i^iEiu,h))+XRNA>f^{E{u,h)). 

In particular, note that 

fi{\E{u,h)\^) = XR«/^i(l-£^(M,h)p)+ XRy^ 2 (|£^(u,h)H, 

P = + Xr"P2, 

m(0) = Xr^Pi(O)+ Xr«P2(0). 

Finally, in order to simplify our notation we set 

GN m,h) := (Gi(«, M, h),..., Gn-i{u, [0], h), G^(u, h))^ 

A{u, h) := {Ai{u, h),..., An{u, h))’^, 

B{u,h) := (BiKh),...,Sw(u,h)f. 


3. The linearized problem 

The above set of equations (2.7) leads to the following associated linear problem 


pdfU — 12 Au + VP 

= / 

in 


divit 

= fd 

in 


-lv{DNUj -1- DjUn)} 

= 93 

on 

K, 

|P] - 2[fiDNUNl - (blXa + crA')h 

/ 

= Pn 

on 


M 

= 0 

on 


dth - UN 

= Ph 

on 

K, 

u\t=0 

= Uo 

in 

R^, 

h\t=o 

= ho 

on 

R^- 


where j = 1,..., - 1 and g = (gi,.. .,gNV- Here, 

P = PiXr^ + P2Xr^ j ^ = ^iXr^ + ^2Xr^ 
with Pi > 0 and Pi > 0 for i = 1, 2. 

The optimal regularity property of the solution of the above problem (3.1) will be of central impor¬ 
tance in the following. To this end, let us recall first the definition of some function spaces. Indeed, let 
m £ N, H C be an open set and X be a Banach space. Then, for 1 < p < oo and s £ R, the Bessel 
potential space of order s is denoted by Hp{il,,X). Moreover, given 1 < p < oo and s £ (0,oo) \ N, 
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Wp{fl,X) denotes the Sobolev-Slobodeckii space equipped with the norm 


|w|(n,x) 




E( 

|a| = W 


||5“/(x)-5' 


7(2/)ir 


A 




|a; _ y|m-|-(s-[s])p 


dxdy) 


i/p 


where [s] is the largest non-negative integer smaller than s and d°‘f{x) = 9l“l/(x)/9x7 ■■■dx^. 
By BUC(p.,X) we denote the Banach space of uniformly continuous and bounded functions on 
and BUC^{Q.,X) denotes the set of functions in C^{Q.,X) such that all derivatives up to order k are 
belonging to BUC{fl,X) for fc e N. For 1 < p < oo, the homogeneous Sobolev space of order 

1 is defined as ijp(ff) = {/ g Lijoc(f^) | ||V/||i^(f 2 ) < oo}. Finally, denotes the dual space 

of Hp,{R^) for 1/p-1- 1/p' = 1. 

The following result due to Priiss and Simonett [22] characterizes the set of data on the right-hand 
sides of (3.1) for one obtains a solution of (3.1) in the maximal regularity space. 

Proposition 3.1. [22, Theorem 5.1], [23, Theorem 3.1]. Let 1 < p < oo, p ^ 3/2,3, a > 0 and 
J = (0, a). Suppose that 


Pi >0, Vi > 0, 7 a > 0 and a > 0, i = 1,2. 


Then, equation (3.1) admits a unique solution {u,0,h) 

U e (ili(J, Lp{R^)) n Lp{J, 

0eLp(J,iji(M'^)), 

[0] e nTp(j, VFpi-i/p(R^-i)), 

h e 1F2 -i/(2p)(j^ Tp(M^-i)) n H^{j, n Lp{J, Wp 3 - 1 /P(M"-^)) 

if and only if the data {f,fd,g,gh,uo,ho) satisfy the following regularity and compatibility conditions: 


f 

e 

Lpij,Lp{mm, 

fd 

G 

H;,ij,Hp\m)) n Lpij,H^pim)), 

g 

G 

(J, Lp{mm n Lp{J, lFpi-i/P(R^-i)))^ 

gh 

G 

ppG-i/( 2 p) Lp{mm) n Lp{J, lTp2-i/P(R^-i)), 

Uo 

G 

Wp2-2/P(R^)^, 

ho 

G 

lFp3-2/P(R^-i), 

div uq 

= 

fm in m, 

[wol 

= 

0 077 i/p > 3/2, 

am 

= 

-Ivi^NUoj + T>jUoJv)l = gjiO) on if p > 3 


for all j = 1,..., N — 1. Moreover, the solution map [(/,/d, p, p/i, mo7o) (^77)] is continuous 
between the corresponding spaces. 
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4. Properties of function spaces involved 


In order to derive estimates for the nonlinear mappings occuring on the right-hand sides of (2.7) we 
study first embedding properties of the functions spaces involved. For a > 0 let J = (0, a) and set 

El (a) = {u€ \ H = 0 }, 

E2{a)=Lp{J,H;,{±^)), 

E3(a) = nLp(J,VFpi-i/P(R"-i)), 

E4(a) = ni7j((J, Wp 2 -i/p(iR^-i)) 

n tFV2-i/(2p)(n Lp{J, VF 3 -i/p(]R^f-i)) 

as well as 

Fi(a) =Lp(J,Lp(R^))^, 

F 2 (a) = n Lp{J,H^0k^)), 

F3(a) = 

F4(a) = {J, Lp{R^-^)) n Lp{J,W^-^^P{R^-^)). 

We then have the following result due to Priiss and Simonett [22]. 


Lemma 4.1. [22, Lemma 6.1]. Let N + 2 < p < oo, a > 0 and J = (0,a). Then the following 
properties hold true. 

a) £ 3 ( 0 ) and £ 4 ( 0 ) are multiplication algebras. 

b) Ei(a) {BUC{J, BUC^fk^)) n BUC{J,BUC{R^)))^ andEi{a) ^ Wp^"^ {J, H^{R^))^. 

c) E 3 {a) ^ BUCiJ,BUCiR^-^)). 

d) £ 4 ( 0 ) BUC\J,BUCk^^~^)) n BUC{J,BUC^{R^-^)). 

e) Wp~^^^^P\j,Lp{R^-^)) r\ H^{J,Wp~^^P{R^-^)) n Lp{J,Wp~^^P{R^-^)) -A £ 4 ( 0 ) 


The crucial point of our proof is the investigatation of the viscosity functions p. To this end, given 
a > 0 , we introduce the function space £ 3 ( 0 ) as 

1 ^ 3 ( 0 ) := {g G BUC{J, BUC{R^~^)) : ||g||jj? 3 (£j) = jiff|| b(7C(j.bc/C(r«”1)) + IsIfsIo) < oo}i 

where | 5 |F 3 (a) = |ff|F 3 (o).i + l 5 lF 3 (a ),2 with 


l5lF3(a),i 




.J Jj 


|t-s| 2 - 


■ dtds 


ijp 


and 


We then obtain the following result. 


\g{t,x') - g{t,y')\P 

\X' - y'\N-2+p 


dx 


'dy'dt^ 


i/p 


Lemma 4.2. Let N + 2 < p < 00, a > 0 , and J = ( 0 ,a). Then the following properties hold true. 

a) £3(0) and £3(0) are multiplication algebras. In addition, 

£3(0) -A BUC{J, BUC{R^-^))^ and £3(0) -A BUC{J, BUC{R^-^)). 

b) If p G BUC^fR) and g G £3(0), then 

ll‘^(ff)llF3(a) ^ I|7 ’I|b( 7C(R) + ll7’||BC/C(R)ls|F3(a)- 

c) There exists a positive constant C such that 

||/ 5 l|F 3 (a) < C'll/l|F3(a) ll 5 llF 3 (a) V / £ ^ 3 {a) and g G £3(0). 




STRONG SOLUTIONS TO TWO-PHASE FREE BOUNDARY PROBLEMS 


9 


Proof, a) The properties for F 3 (a) is essentially given in Lemma 4.1 (a) and (c). The embedding 
^ 3 ( 0 ) ^ BUC{J,BUC{'M.^~^)) follows from the definition of F 3 (a). We thus only show that F 3 (a) is 
a multiplication algebra. For f,g G F 3 (a) it follows that 

II/S'I|bUC(J.SC/C(R'^-i)) ^ II/I|bC/C(J,B(7C(R'V-i))||S'||b(7C(J.BUC(R'^-i)) ^ ll/llF3(a) 11511^3(0)’ 

Considering | • |F 3 (a),i we see that 

. /■ /■ ... M/P 

lffflF 3 (a),l^llfllBUC(J,BUC(R^-^))(/ / ---TITT- dtds)' 

J.I JJ R — 5 2+7 


+ I|5||bC/C(J,B(7C(R'V-1)) 


\t-. 

II/W-/(s)IIl,(R^-1) . , \1/p 

dtds\ 


IjJj |t-s|2 + 2 

- ll/llF3(a)l5lF3(a),l + llff llF 3 (a) I / IPsCa),!' 

Similarly, |/g|F 3 (a). 2 ^< I|/Ili 3 (a) | 5 |F 3 (a ).2 + || 3 |lF 3 (a) l/lF 3 (a). 2 - This yields ||/ff||F 3 („) < C'll/llF 3 (a) ll 5 llF 3 (a)> 
which implies that F 3 (a) is a multiplication algebra, 
b) By the mean value theorem 

lArfIp.,., = ( J J --TFI-<«*) 


+ 


\ip{g{t,x')) - (p{g{t,y'))\P 


JJR«‘-idRN-l \x'— y'\^ 


\ J 

’dy'dt] 




+ 


it- 

\g{t,x’) - g{t,y')\P 


J jR^v-i \x' — y'\^ 


wdty^} 


< II'/5||b[/c<(r)|5|f3(q 


which yields the required inequality, 
c) Obviously, 

ll/5llLp(J.Lp(R«‘-i)) < II/IILp(J.Lp(R'^-i))I|5||b(7C(J,BC/C(R'^-i)) < ll/llF3(a) ll5llF3(a)- 

On the other hand, we see that by Lemma 4.2 (1) and by calculations similar to a) there exists a 
constant C > 0 such that for i = 1, 2 

l/5lF3(a),i < ll/llBUC(J,B(7C(R'V-i))|3|F3(a).i + IIsIIb( 7C( J,BUC(R«‘-i)) I/IfsIoIR - C'H/llFs/a) lls llF3(a) ’ 
which combined with the above inequality completes the proof. □ 


We next recall basic properties of functions which are Frechet differentiable. Let X and Y be 
Banach spaces and U C X he open. We then denote the Frechet derivative of a differentiable mapping 
$ : {7 —/ y by : U -G C(X,Y) and its evaluation for it G 17 and v G X hy [7?$(u)]ii G Y. 
Moreover, a mapping $ : 17 —>■ 1" is called continuously Frechet differentiable if and only if $ is Frechet 
differentiable on U and its Frechet derivative is continuous on U. The set of such continuously 
Frechet differentiable mappings from 17 to F is denoted by C^{U,Y). 

In the sequel, we will make use of the chain and product rule for Frechet differentiable functions. 
In fact, in addition let Z be a further Banach space and suppose that the mappings / : 17 —/ F and 
g : Y ^ Z are continuously Frechet differentiable. Then the composition F = gof: U^Z is also 
continuously Frechet differentiable and its evaluation at x G 17 and x G 7f is given by 

[DF{x)]x = [Dg{f{x))][Df{x)]x. 
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For the product rule, suppose that there exists a constant M > 0 such that for every y €Y and z G Z 

\\yz\\Y < Mllyllrllzllz, 

and also that f : U and g :U ^ Z are continuously Frechet differentiable. Set F{x) = f{x)g{x) 
for X G U. Then F : U ^ Y is also continuously Frechet differentiable and its evaluation a.t x G U and 
a: G X is given by 

[DF{x)]x = g(x)[Df{x)]x + f{x)[Dg{x)]x. 

Now, we define the solution space E(a) and the data space F(a) for a > 0 by 

E(a) := {(u, 9, tt, h) G Ei(a) x E 2 (a) x E 3 (a) x E 4 (a) | |0J = tt}, 

F(a) := Ei(a) x F 2 (a) x F 3 (a) x F 4 (a). 

The spaces E(a) and F(a) are endowed with their natural norms, i.e. 

||(U, 0,7r,/l)||]E(a) : = ||u||Ei(a) + ll^'l|E2(a) + lkl|E3(a) + ll^llEifa), 

\\if,fd,g,gh)\\¥{a) ■ = ||/llFi(a) + ll/d||F2(a) + ll5llF3(a) + lISli IIf 4 (o) ■ 

Finally, we consider for (u, 9, tt, h) G E(a) the nonlinear mapping N which is defined as 

(4.1) N{u, 9, TT, h) := (F(m, 9, h),Fd{u, h),G{u, tt, h),Gh{u, h)), 

where the terms on the right hand side are defined as in Section 2. For functions u = (ui,... 
dehned on we set 

(4.2) : = (ul,... Uj:=XR^Uj, 

u : = (Ui,... ,M^), 

Recalling the definition of F(u,h) and £(u,h) in (2.6), the following lemma and its corollary shows 
that various functions occuring in the definition of N in (4.1) are Frechet dfferentiable. 

Lemma 4.3. Let IV + 2 < p < oo, a > 0 and J = (0,a). Then the following assertions hold true. 

a) For (j) G BUC^{M.), the mapping 

(fi : BUG{J,BUG{]k^)) BUC{J, BUCiR^)) 

is continuously Freehet differentiable. 

b) For ip G BUC^fR), the mapping 

if : F 3 (a) ->• F 3 (a) 

is continuously Frechet differentiable. 

c) Let (p G BUG^(M.) and for u defined as in (4.2) set 

<^‘‘{u,h) := (p{\E{u‘‘,h)\'^) ford =1,2. 

Then : Ei(a) x E 4 (a) —> BUG {J, BUG {M.^)) is continuously Frechet differentiable. 

d) Let Ip G BUC^{M.) and for u defined as in (4.2) set 

4''^(u, h) := 'ip{\'yoE{u‘^, /i)p) for d=l,2, 

where 70 denotes the trace to Rq . Then 4''* : Ei(a) x E 4 (a) —>■ F 3 (a) is continuously Frechet 
differentiable. 

Proof, a) We show first that the mapping p is Frechet differentiable. To this end, let /, / G Z := 
BUG{J,BUGi^^)). Then 

tU + /) - Pif) - T{f)f = [ {pif + &f) - tU)) d9f, 

Jo 
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which implies 

yif+f) -^if)-v{f)f\\z/\\f\\z < f yif+of) - (pif)\\zd0. 

Jo 

Since ip S BUC{B.), the term on the right hand side above tends to 0 as WfllsuciJBUC{m’^)) 

Thus [Dip{f)]f = <p{f)f. Next, we show the continuity of the Frechet derivative at /o G Z. For h € Z 
we have 

\\D(p{fo + h) - D(p{fo)\\c(z) = sup \\[D(p{fo + h)]f-[Dip{fo)]f\\z 

ll/llz=l 

= sup \\p{fo + h)f -p{fo)f\\z 

\\f\\z = l 

< W^Pifo + h) -p{fo)\\z, 
which tends to 0 as \\h\\z —>■ 0 since p G BUC{^). 

b) For /, / G F 3 (a) we obtain 

^(/ + /) - ^(/) - ^(/)/ = / V - + 0f)ffde. 

Jo 

By Lemma 4.2 a) and b) 

ll'0(/ + ^/)/lli3(a) <C'{||^I|bC/C(R) + II^'||bC/C(R)I/ + ^/|F3(a)}ll/lli3(a) 

<C’(l+ll/llF3(a)+^ll/li3(a))ll/l|3(a)^ 

which implies that [Dtjj{f)]f = tp{f)f. Next, in order to show the continuity of the Frechet derivative 
at /o G F 3 (a), let h G F 3 (a). Then 

+ = sup^ ll[-Di/’(/o + M]/-[^V'(/o)]/||f3(a) 

Il/lb3(a)=l 

<cmfo + h)-p{fo)\\i^^^^ 

by Lemma 4.2 a). Since ip G BUC'^pR.), Lemma 4.2 b) implies that ip{fo + h) and ip{fo) are in F 3 (a). 
Taylor’s formula and Lemma 4.2 a) yield 

mio + h)- l/-(/o)|lr3(,) ^ ^ + ^^)llF3(a) ^^ll^llF3(a)- 

'J 0 

The latter terms tends to 0 as ^ 0, which completes the proof. 

c) By Lemma 4.1 b) and d), the mappings 

(4.3) (u, h) i-A E'^{u, h) : Ei(a) x E 4 (a) ^ a,s a; i-A \x\^ : —>■ Z 

are continuously Frechet differentiable. The chain rule thus yields that 

El (a) X E 4 (a) Z : {u, h) \E{u‘^, /i)p 

is continuously Frechet differentiable, too. Applying assertion a) and the chain rule again implies that 
: El (a) X E 4 (a) Z is continuously Frechet differentiable. 

d) Note that for (u, h) G Ei(a) x E 4 (a) and i, j = 1,..., A^, we obtain 

(4.4) \\{jQdiu'p,dih,didjh)\\f^f^^^ + \\{-fodiu‘^,Vh,diVh)\\w^(a) < C{a,p)\\z\\^(^a) 

In fact. Lemma 4.2 b) and [18, Theorem 4.5] for s = 1/2, m = 1, and /r = 1 yield 

+ ||7o5iU‘^||F3(a) <ll7o3iU^||B(7C(J,BUC(R'^-i)) + ll7o9iM‘^||F3(a) 

<C(||a*U^||uClC'(J.B!7C(Ri/)) + ll^»“'*llw/''"(J,Lp(R«))nLp(J.Hi(R«))) 
^C'll'ullE^(a), 
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which implies the required properties of u. Concerning h, the desired properties follow from the 
definition of E 4 (a) and Lemmas 4.1 and d). By (4.4) and Lemma 4.2 a) the mappings 

{u, h) n- 'yoE{u‘^, h) : Ei(a) x E 4 (a) —>• F 3 (a)'^^'^ and x i-A \x\^ : F 3 (a)'^^^ ^ F 3 (a) 

are continuously Frechet differentiable. By the chain rule 

(u, h) i-A \^qE{u'^, h)\^ : El (a) x £ 4 ( 0 ) —>■ £ 3 ( 0 ) 

is continuously Frechet differentiable. Finally, combining this with assertion b), the chain rule yields 
the assertion. □ 

For $ S BU C'^(K), u as in (4.2) and i, k,l,m,q,r = 1,..., and d = 1, 2 we now introduce the 
functions 

:= ^{\E{u^h)\^)E,,{u^h)Euiu‘^,h)drnd,4, 

^ijkipqriu^h) := ip{\E{u‘^,h)\'^)E,j{u‘^,h)Eki{u‘^,h)Emqih)uf, 

:= {ip{\E{u‘^,h)\^) - ip{0))didjui, 

:= i‘P{\E{u‘^,h)\^) - ip{0))Eij{h)ut, 

as well as for ip € BUC^{M.), u as in (4.2), i,j,k = l,...,N and d = 1, 2 the functions 

'^iiu,h) := {'0(|7o-E(m‘^,/i)P) -'0(O)}7o9iu'^, 

■■= tp{\-foE{u‘^,h)\‘^){-fod,u‘^)djhdkh. 

Corollary 4.4. Let N + 2 < p < 00, a > 0 and J = (0,a). Then the following assertions hold true. 

a) Assume that i,j.,k,l,m,q,r = 1,...,A^ and d = 1,2. For (j) £ BUC^iM.) and u as in (4.2) the 
mappings 

Afifc ■ ^1 (o) X £4(0) Lp{J, LpfR^)) 

are continuously Frechet differentiable. Moreover, their values and their Frechet derivatives at 
{u, h) = ( 0 , 0 ) vanish. 

b) Let i,j,k = 1,..., N and d = 1,2. For tp £ BUC^{M.) and u as in (4.2) the functions 

vE-f,0f^.,Sf^.fc:Ei(a)xE4(a)^F3(a) 

are continuously Frecht differentiable. Moreover, their values and Frechet derivatives at (u, h) = 
(0,0) vanish. 

Proof, a) We only prove the assertion for ^fjkimqr'^ remaining assertions can be proved in an 
analogously. By (4.3), Lemma 4.3 a) and the product rule, the function 

{u,h) HA ip{\E{u‘^,h)\^)E„{u'^,h)Eki{u‘^,h) : Ei(a) x £ 4 ( 0 ) ^ BUC{J, BUC{t^)) 

is continuously Frechet differentiable. Moreover, 

u HA dmdquf : £i(a) Lp{J, Lpfk^)) 

is continuously Frechet differentiable, which combined with the above assertion and the product rule 
applied to the situation X = Ei(a) x £ 4 ( 0 ), Y = Lp{J, Lp(R^)) and Z = BUC{J, BUC{R^)) implies 
that ^fjkimqr ' ®'i(®) ^ ]Ei(a) Lp(J, Lp(]R'^)) is continuously Frechet differentiable. In addition, it 
is clear that 0 ) = 0 and ( 0 , 0 ) = 0 . 

b) We only prove the assertion for Since 

(4.5) u HA 'yodiu'^ : £i(a) £ 3 ( 0 ) or £ 3 ( 0 ) 
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is continuously Frechet differentiable, it follows from (4.4) that 

u i-A '!/'( 0 ) 7 o 9 iu'^ : El (a) —F 3 (a) 

is continuously Frechet differentiable, too. On the other hand by Lemma 4.3 b), (4.5) and the product 
rule applied to X = Ei(a) x £ 4 ( 0 ), Y = £ 3 ( 0 ) and Z = £ 3 ( 0 ), the mapping 

(u, h) i-A i/;(| 7 o£'(u'^, h)|^) 7 o 9 iu‘^ : Ei(a) x £ 4 ( 0 ) £ 3 ( 0 ) 

is continuously Frechet differentiable. Finally, 4'f(0,0) = 0 and thanks to product rule 

D^f{u,h) = -fod^u‘^[D{tlj{\'yoE{u‘^,h)\'^) - •0(0))] + {tlji\'yoE{u‘^,h)\‘^) - il;{0))[D{-fodzu‘^)], 

which implies that £)4'(*(0,0) = 0. The proof is complete . □ 


5. The nonlinear problem 

Let us recall from (4.1) that for {u,6,tt, h) G £(a) the nonlinear mapping X was defined as 
N(u, 0, TT, h) = (E(u, 9, h),Fd{u, h),G{u, tt, h),Gh{u, h)). 

We start this section by examining properties of the nonlinear mapping N. 

Lemma 5.1. Let N + 2 < p < 00 , a > 0 and r > 0. Suppose that Pd{s) G ^^([O, 00 )) for d = 1, 2 and 
in addition that pi > 0, p 2 > 0, 70 > 0 and tr > 0 are constants. Then 

N G C'^(BE(a)(T),£(a)), iV(0) = 0 and DN{0) = 0. 

Proof. We treat here in detail only the terms A{u, h) and i3(M, h) which are defined in Section 2. The 
remaining terms may be treated as in [22, Proposition 6.2] and [23, Proposition 4.1]. 

The term A(u, h): 

Let (u, 9, TT, h) G i?E(a)(^) and recall that Ai is given for * = 1,..., by 

N 

A,{u, h) = '^ {Af^’\E{u, h)) - Al’^’\0)){djdkui + djdiUk) 

N 

+ ^ {Ai’^’\E{u, h)) - Ai’>^’\0)){E,,{h)ui + Ep{h)uk), 

where Ejk{h) are defined as in (2.5) and as in Section 2 by 

Af^’\E{u,h)) = XRi;^A{f\E{u,h))+XR^A^f^/{E{u,h)), 

A^f^/{E{u, h)) = i (2fidi\E{u, h)\‘^)Eij{u, h)Eki{u, h) + pd{\E{u, h)\‘^)Sik6ji^ 

for d = 1, 2. Observe that Af^’\E{u, h)) may be represented as 

.. 2 

Af^'^Eiu, h)) = -J2 ( 2 Ad(|-B(w^ h)\^)E,j{u‘^, h)Eki {u\ h) + pd{\E{u‘^, h)\^)5,uSji) , 
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and thus 


2 N 

Mu,h)=Y, J 2 + 

d—1 —1 

2 N 

d—1 j,k,l—l 
2 N 

+E E h)\^)E,,{u‘‘, h)Eki{u^, h){Ejk{h)uf + E,i(h)ui) 

d—1 _7,fc,/—1 
2 N 

d—1 _7,fc,/—1 

for i = 1,..., By Corollary 4.4 a) 

^ e C'^(SE(a)(r),Fi(a)), y 4 ( 0 , 0 ) = 0 and £>^( 0 , 0 ) = 0 . 

The term B{u, h): 

Let (m, 9 , TT, h) G i 3 E(a) (''’)■ By Lemma 4.1 b) and d), each term appearing in B{u, h) is continuous with 
respect to the space variable. In particular, this implies 

7o/id(|£(u‘^, h)p)) = (j.d{\loE{E^, h)\'^)), lo{{dNU%){djh)} = {'jodNU%){-jQdjh). 

Thus, B{u^ h) may be rewritten as 

2 

Bj{u,h) = - fid(.\7oE{u'^, h)\'^){jodNU%)djh 

d^l 

2 

+ '^(-^fiBd{\loE{u‘^, h)p) - fid{0)){lodNUj + 7odjU%) 
d^l 

2 N-1 

-EE {-^ytJ-d{\'loE{u‘^,h)\'^){jodjui + -fodkUj)dkh 

d—1 k—1 
2 N-1 

+EE {-Nnd{\loE{u'^,h)\'^){dkhjodNu'^ + djh^QdNui)dkh, 

d—1 k—1 
2 

BN{u,h) =2Y,i-^f{Td{\loE{u‘‘,h)f) - ^idmiodNui 
d^l 

2 N-1 

+EE {-Nnd{\loE{u'^,h)\'^){dkh)'^^QdNU% 

d—1 k—1 
2 N-1 

-EE {-Nfid{\loE{u‘^, h)p)(7o5ArMfc + 'yodkU%)dkh. 

d^l k^l 

These representation combined with Corollary 4.4 b) yields 

B G C\BEia)ir),¥3ia)), S( 0 , 0 ) = 0 and ££( 0 , 0 ) = 0 , 
which implies the assertion. 

Finally, we return to the nonlinear problem ( 2 . 7 ). We define the space of initial data I by 

2-2/p f!nN\N „ XX/3-2/pfTnN-l 


□ 
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and define also for z £ E(a) and (ito, /lo) € I the mapping $ by 

$(z) := L~^{N{z),uo,ho). 

Here L is defined by the left-hand side of the linear problem (3.1) with u = /r(0). Observe that the 
invertiblility of L is guaranteed by Proposition 3.1 since > 0 for z = 1,2 by assumption and 

N{z) € F(a) for z G E(a) by Lemma 5.1. The following result shows that the problem (2.7) on the 
fixed domain admits a unique strong solution provided the data uq and ho are sufficiently small in 
their corresponding norms. 

Proposition 5.2. Let N + 2 < p < oo and a > 0. Suppose that pi € C'^([0, oo)) for i = 1,2 and that 

Pi > 0, Pi{0) >0, 7 o > 0 and cr > 0. 

Then there exist positive constants Eq and So (depending on a and p), such that system (2.7) admits a 
unique solution (u,9,h) in B^i^^p^^So) provided that the initial data {uo,ho) £ I satisfy the compatibility 
conditions 

(5.1) Ip{\E{uo, /io)|^)£'(uo, ho)no - {uq • p{\E{uo, ho)\‘^)E{uo, ho)no}no] = 0 on 

lOoj = {no ■ p{\E{uo, /io)P)L;(uo, ho)nol + {[ph + aA')ho - aTLiho) on , 

divito = Fd{uo, ho) in R^, |uo| =0 on Rq 
as well as the smallness condition ||(uo,ho)||i < Eq- 
Remark 5.3. The compatibility conditions (5.1) are equivalent to 

-lp{0){dNUj + djUN)l = Gj{uo,l0oj,ho) on 
| 6 >o] - 2|/Lt(0)5jvuoJvl - (Mba + crA')ho = GAr(uo, ho) on , 
divuo = Ed{uo, ho) in |itol = 0 onR^ 

for j = 1,..., iV — 1 and G = (Gi,..., Gn)^ defined as in Section 2. 


Proof. Observe that DN is continuous and DN{0) = 0 by Lemma 5.1. Thus, we may choose Jq > 0 
small enough such that 


sup < 

ZGB]s,(^a) (25o) 


1 

2||T“^||z;(F(a)xI,E(a))’ 


where r > 0 is a sufficiently large number. For z £ B^(^a)iSo), the mean value theorem implies 


||$(^)||E(a) < ||i"^||c(F(a)xI.E(a)){||-/V(z)||F(o) -I- ||(mo,^o)||i} 

= \\B ^||c(F(a)xI.E(a)){l|^(2:) - iV(0)||f(a) -|- Eq} 

<\\B ^||c(F(a)xI.E(a)){(_ SUp ||-DiV(z) || (r) .F(o)) ) II ^ |lE(a) + } 

2GBE(a)(<5o) 



£o\\L ||£(F(a)xI,E(a))- 


Choosing Eq in such way that 0 < Eq < <5o/(2||L ^ ||£(F(a)xi.E(o))), we obtain ||$(z)||E(a) < <^o- Hence, 
$ is a mapping from i?E(o)(^o) into itself. 

Let zi,Z 2 £ i?E(a)(^o)- Noting that $(z:i) — d)(z 2 ) = L“^(A^(zi) — iV(z 2 ), 0 , 0 ), we obtain by the 
mean value theorem 


||$(zi) - ^'(2:2)||E(a) < ||.^ ^|U(F(a)xI.E(a))||A^( 2 :i) “ ^(z2)||f(o) 

< l|-^”^llc(F(a)xI,E(a))( SUp || £>^^( 2 ) || (r) .F(a)) ) ) 1^1 “ ^^2 ||E(a) 

(2(5o) 

< 2 ~ ^ 2 ||E(a), 
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which implies that <i> is a contraction mapping on i3E(a)(<^o)- The contraction principle yields the 
existence of a unique solution of (2.7) in 73^(0) (<^o)- D 

Proof of Theorem 1.1. Note that the compatibility conditions of Theorem 1.1 are satisfied if and only 
if (5.1) is satisfied. The mapping &ho given by 

■= + hoiO) for G 

defines for hg G Wp ^^^(R^“^) a C^-diffeomorphism from onto fl(0) with inverse ©/(J^(a;', xjv) = 
(x',XN — ho{x')). Thus there exists a constant C{ho) such that 

C{ho) < (7 (/IQ ) 11 UQ 11 jy 2 - 2 /p ^ ^ . 

Hence, the smallness condition in Theorem 1.1 implies the smallness condition in Proposition 5.2. 
Proposition 5.2 then yields a unique solution {u,9,h) G BE(a){So) of (2.7). Finally, setting 

(v, Tt) = (0*M, 0 *7r) = (u O 0“^, TT O 0“^), 

where 0* is defined as in (2.4), we obtain a unique solution (y, tt, h) of the original problem (1.2) with 
the regularities stated in Theorem 1.1. The proof is complete. 

□ 
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